Quantum scattering theory on graphs with tails 
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We consider quantum wallcs on a finite graphs to which infinite tails are attached. We explore 
how the propagating and bound states depend on the structure of the finite graph. The S-matrix 
for such graphs is defined. Its unitarity is proved as well as some other of its properties such as its 
transformation under time reversal. A spectral decomposition of the identity for the Hamiltonian 
of the graph is derived using its eigenvectors. We derive formulas for the S-matrix of a graph under 
certain operation such as cutting a tail, attaching a tail or connecting two tails to form an edge. 
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I. INTRODUCTION 

Quantum walks have become important tools for mod- 
eling and analyzing the behavior of a quantum system in 
quantum information theory. Two seemingly different 
types of quantum walks have been defined — discrete- 
time and continuous-time quantum walks. Discrete-time 
quantum walks come in several flavors — on a regular 
undirected graph equipped with a coin space in on 
undirectedgraphs in [3, y] , on the edges of an undirected 
graph in H, |5[ or corresponding to a classical Markov 
chain in [6 . Continuous-time quantum walks first ap- 
peared in 7|. Despite the similarities in their behav- 
ior, a way to obtain one as a limit of the other, some- 
thing easily done in the classical case, was not known. 
This problem was addressed in P. It was also estab- 
lished that continuous-time quantum walks on sparse, 
low-degree graph are universal for quantum computation 
i- 

Quantum walks, similarly to random walks in the clas- 
sical case, are very useful for developing quantum algo- 
rithms. Some of the first algorithms discovered employ- 
ing this approach are element distinctne ss |lOl| . matrix 



product verification [11 
commutativity testing 



triangle finding [12| | and group 
13| . The quantum algorithm for 
the glued-tree graph is based on a continuous-time quan- 
tum walk and was proven to be exponentially faster than 
its classical counterpart in [l^ . These examples of quan- 
tum algorithms are based on the fact that the quantum 
walks "hits" a special vertex polynomially or exponen- 
tially faster when compared with the classical walk or 
the fastest known classical algorithm for the problem in 
question. Thus the question for a proper definition of a 
"hitting time" for quantum walks arises. In the discrete- 
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time case definition for the hitting time were given and 
explored in [l^ [l^, 113] , and for the continuous-time case 
in [isj . An optimal quantum algorithm for evaluating a 
balanced NAND tree proposed in [l^ is another exam- 
ple of a continuous-time quantum algorithm for which 
the value of the NAND tree depends on whether a cer- 
tain part of the graph is hit or not. The premise for 
the graph is somewhat different though — it involves the 
consideration of a graph with infinite tails attached to 
it. This result spurred a plethora of quantum algorithms 
for evaluating formulas |28l . [29j . Scattering theory for 
discrete-time quantum walks was developed in [l^, HH ■ 
Another type of quantum walk-based algorithm involves 
generating a saniple from the uniform distribution over 
a graph [l^, |2^, [23| . Thus one needs to define "mixing 
time" for quantum walks. Properties of mixing times and 
lower bounds are proven in [25, 26, 27]. 

In this paper we study the model considered in p^ . 
There, two tails of infinite length are connected to a finite 
graph that contains information about the values of the 
variables that enter the NAND formula. A wave with 
fixed energy (or at least with narrowly peaked energy 
spectrum) is sent along a tail. After a time T equal 
to square root of the number of variable N the walk is 
measured to be found on one of the two tails as the wave 
is either totally reflected or transmitted for the selected 
energy. As there is a direct correspondence between the 
value of the reflection coefficient and the value of the 
NAND formula we deduce the value of the formula. 

As we already noted, |9j] proves that continuous-time 
quantum walks on graphs with infinite tails are univer- 
sal for quantum computation. In this model, as in the 
previous one, the energy of the waves incoming along 
the tails is fixed. An arbitrary state of a qubit is rep- 
resented by superposition of waves incoming along two 
tails. Thus N qubits are represented by 2^ tails. A one- 
qubit (two-qubit) gate is a graph that has 2 (4) incoming 
and 2 (4) outgoing tails. Universality follows from the ex- 
istance of graphs that implement the control- not, phase 



2 



and Hadamard gates. 

The above two papers and numerous other algorithms 
based on this model prompt us to explore scattering the- 
ory on graphs in continuous time. In section |lTl we give 
the basic model for the graph with infinite tails on which 
the quantum walk evolves and describe both propagat- 
ing and the two kinds of bound states for the form of 
Hamiltonian we consider. In section Hill we prove the or- 
thogonality between propagating and the bound states. 
The S-matrix is defined and its unitarity proven in sec- 
tion |TVl In the following section |V] we explore the rela- 
tionship between energy eigenstates and vertex states. In 
sections [VIIVIII and I Vllll formulas for the S-matrix under 
the operations of cutting, attaching and connecting tails 
are derived. In the last section IIXI the unitarity of the 
S-matrix under these operations is proven. 



II. ENERGY EIGENSTATES 

We want to introduce scattering on an undirected 
graph in continuous time. With each vertex w of a graph 
G we associate a normalized state \v) in a Hilbert space 
Tic, the graph Hilbert space, with the property that 
states corresponding to different vertices are orthogonal 
to each other: {vi\vj) = Sij. The structure of the graph 
determines the unitary evolution in the graph Hilbert 
space. The Hamiltonian we consider is a Hermitian op- 
erator on the graph Hilbert space which in the basis of 
vertex states is given by minus the adjacency matrix of 
the graph [H]. 

In order to investigate the scattering properties of 
graphs we need to connect "tails" to a finite graph. A 
tail is a semi-infinite linear graph with its end connected 
to one vertex of the graph. More than one tail can be 
attached to any one vertex of the graph. The vertices of 
each tail are number from 1 to infinity and the vertex to 
which it is attached is labeled with a 0. To identify dif- 
ferent tails sometimes we will label them with the vertex 
of the finite graph to which they are attached. If more 
than one tail is attached to the same vertex we can use a 
second index to differentiate between different tails. The 
Hamiltonian on a tail is defined as in the above as minus 
the adjacency matrix of the tail. Explicitly 

oo 

iJ*«' = -^(|n)(n+l| + |n + l)H). (1) 

n=l 

Let's consider a linear graph infinite in both directions. 
It's easy to see that the energy eigenstates are given 
by waves propagating in the positive or negative direc- 
tion: {n\k)± — e^**^" with k S (0,7r). States that cor- 
respond to fc = or fc = TT do not behave exactly like 
propagating states and need to be analyzed separately. 
The energy corresponding to a state with absolute value 
of its momentum given by |fc| regardless of propagation 
direction is E±k — — 2cosfc. 



A. Propagating states 

We want to find the continuous energy spectrum and 
propagating energy states of a graph with tails. We start 
by considering the usual stationary Schrodinger equation: 

H^\^)^E\i;). (2) 

The Hamiltonian is given by 

H^=Hg+J2Y. (^^™ - |0.)(l.r„| - |l.rn)(0.|) . 

(3) 

In the above equation the first sum is over all vertices 
of Hg to which tails are connected and the second sum 
is over all tails connected to a certain vertex v. The 
number of tails connected to vertex v is denoted by my 
and total number of tails by m = . The notation 

\nvm) refers to the n-th vertex of the m-th tail connected 
to the vertex v of the graph G and as noted above |0„) 
is the state associated to the vertex v, in other words 
|0.) = Iz;). 

Let's denote by \k, vm), with k G (0, tt), the propagat- 
ing energy eigenstate of the above Hamiltonian which has 
an incoming component only on the m-th tail connected 
to the vertex v of the graph G. We denote the restriction 
of this state on the graph G by |fc,i;m)'^. Then 

{nym\k, vm) = e-'''^" + w(fc)e*". (4) 

This state on all of the rest of the tails will look like: 

(nv'm'\k, vm) — ty'm' ,vmik)e . (5) 

Explicitly 

oo 

\k,vm) = \k,vmf + ^(e-^''^" -f w(fc)e''^") |n„„) 

n=l 

I ^ 

(v',m') n=l 

The prime on the sum means that we sum over all or- 
dered pairs (w',m') different from (w,m), a convention 
we will use further in the paper. The state \k, vm) de- 
notes the restriction of |fc, vm) on the finite graph G. We 
note that these states are not normalizable as propagat- 
ing states should be. Substituting the above expression 
for them in ^ and taking into account that the energy 
corresponding to this state is Ek = —2 cos k, after some 
cancelations we obtain: 

Hg \k,vm)^ + (1 + r„,„) - {e-''' + ry^e''') \0v) 

— ^ {Oy\k,vm)'^ l^v'm') — ~-2cosk \k,vm)^ . 

(d' ,m') 

(7) 
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As the terms that hve on the tails should cancel it follows 
that 

rvm + 1 = {Qv\k, vm)'^, 

tv'm',vrn = {Ov' \k , Vm)'^ . (8) 

Substituting these back in equation ([7]) it becomes 

Hg \k,vmf - (e-^'^ - e^'^) |0„) - e^'= |0„) (0,|fc)^„ 
-e''' \0'^) {0'Jk,vm)^ ^ -2 cos k\k,vmf . (9) 

{v' ,m') 

Finally the equation that \k^vm) satisfies is 



+ 2 cos - e'^ \^v'){^v' I ] \k, 

\ v' / 



= (e-*'^- - e''^) |0.) . (10) 

Another, more compact form of the equation that will 
turn out to be useful, we can get if we use the variable 
z = e**^. Then 

(/ + zHg + z^Q) \z, vmf = (1 - z2) |o„) , (11) 

where we introduce the following operators 



i? = ^rrv 1 0,.}(0, 

v' 

Q = I -R. 



(12) 
(13) 



For convenience we denote A{z) = I + zHc4 + z^Q. Later 
we will prove that when z lies on the unit circle in the 
complex plane equation (jlip will always have a nonzero 
solution, and thus a propagating state will exist with 
reflection and transmission coefficients defined by ([5]). 
Here we want to note that when k € (— tt, 0) the states 
\k,vm) = \ — \k\,vm) are defined exactly as above. They 
can be expressed as linear combinations of the states 
\\k\,vm). The formula for this will be given in Chapter 
IVl Another thing we want to note is that the operator 
A{z) and the solution to pT|) associated with it is an 
instance of the so-called quadratic eigenvalue problem. 
For a good overview of the subject with many applica- 
tions and examples discussed one can look at [30]. The 
equation for the bound states derived in the next chapter 
is exactly the equation for the the right eigenvectors of 
A{z). 

The two cases k = and k = tt need to be analyzed 
separately. The states corresponding to those values of k 
still need to satisfy equation (jlip for z = 1 and z = — 1, 
respectively. Thus 



{I±Hg 



Q)i±i)!L 



= 0. 



(14) 



This equation may and may not have solutions, which is 
in contrast with the other propagating states which al- 
ways exist as we shall prove below. If a solution |e}^„ 



with e = ±1 exists and Q le)^^ 7^ then such a state 
will be nonzero on some of the tails and thus it will not 
be square-summable. A difference from the rest of the 
propagating states is that there may not be linearly in- 
dependent solutions for each tail. The notion of incoming 
and outgoing waves is also lost when it comes to these 
states. Thus refiection and transmission coefficients can- 
not be defined properly. In some respects these states 
behave more like bound states in the sense that they sat- 
isfy an equation identical to the bound state equation 

m- 

In the case when Q |e)„„j — the full state |e}^„i will 
be zero on each tail and thus be a real bound state of the 
second kind defined below. 



B. Bound states 

Bound states are solutions to ^ which are normaliz- 
able. This is possible if and only if the amplitudes of such 
bound state on the tails of the graph are either exponen- 
tially decaying or zero. Thus two kinds of bound states 
exist. The first are states for which there is at least one 
tail on which the state is non-zero. The second kind of 
bound states are zero on all tails of the graph. We will 
see that the energy of these two kinds of bound states is 
qualitatively different. 



1. Bound states of the first kind 

Let us denote by \>Cf,) a solution of ([2]) normalized such 
that on the tails of the graph it has the form 



\^b) 



(15) 



As we want this to be a bound state of the first kind, 
ctvm 7^ for at least one pair {v, m). We have the follow- 
ing expression for the bound state 



(16) 



71=1 



with N^^ being a normalization factor. Here again l^Cf,) 
is the part of the state living on the finite graph G. By 
applying the tail Hamiltonian ([T]) to the tail of this state 
for which a„m ^ we find its energy E'^^ = — (e'^'' -I- 
g-x6^_ poi- the energy to be a real number e^^*" should 
be either a real number or a complex number that belongs 
to the unit circle. From the requirement that the state 
is normalized 



(17) 
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it follows that the infinite sum J2^=i 1^ is conver- 

gent which leads to the condition that |e^^''|^ < 1. This 
means that for bound states of the first kind — e^^' 
should be a real number with absolute value strictly less 
than 1 or in other words ^{xh) > and = or tt. 

Inserting \>ct) in ^ leads to the following equations: 

{Hg + 2cosh>f6 - e-'^'R) \xbf = 0, 

If we again use a change of variables, = e^'^'', the 
equations take the following form: 

{l + ZbHG + z^Q)\zbf ^0, (18) 
avm{zb) = (0„„|zfc)'^. (19) 

We note that the operator A{z) defining the propagating 
states through (jlip appears in the definition of the bound 
state of the first kind in the above equation but with 
z = Zb real and satisfying |z6p < 1. 

To determine the normalization factor A^^^ we consider 
equation (|17p in the new variable taking into account 
PU)) . After some simplifications we obtain 

\N,,\^ = ^-^^ ^. (20) 

<^{zb\I-zlQ\zbf 

2. Bound states of the second kind 

In contrast with the bound states of the first kind the 
bound states of the second kind have zero overlap with 
any vertex on any tail. They live entirely on the graph 
G. From this it easily follows that such a bound state \(3) 
has the property {v\(3) — for any vertex u £ G to which 
a tail is attached. Because of that equation ([2]) for a such 
bound state reduces to 

Hg 1/3) = Eb \P) . (21) 

We see that such a state is an energy eigenstate of the 
finite graph G with energy Efj . We will show in the next 
section that if such bound states exist and their energy 
is less than 2 this will lead to non-invertibility of the 
operator I+zHg+z^Q for some z on the unit circle. This 
will necessitate the redefinition of propagating states for 
the value of z. 



III. ORTHOGONALITY OF PROPAGATING 
AND BOUND STATES 

In this section we will prove the existence of propagat- 
ing states for any k £ (0, tt). The bound states of the 
first kind are obviously orthogonal to any propagating 
state because their energy, E^^ — — 2cosh>«:f,, is always 
greater in absolute value than 2 and the energy of the 
propagating states, Ek = — 2cosfc, is always less than or 



equal in absolute value to 2 (it is well known that eigen- 
states of a Hermitian operator with different eigenvalues 
are orthogonal to each other). 

Let's assume now that the operator A{z) = I + zHg + 
z^Q is not invertible for some z = zq on the unit circle. 
Then there is a normalized state \u) living totally on the 
finite graph such that 

(/ + z^Hg + zlQ) \u) = 0. (22) 

Multiplying the above equation on the left with {u\ and 
denoting {u\Hg\u) = G M and {u\Q\u) = 1 — 
{u |(X;„ |0i,)(0i,|)| m) = 1 - r e M we ge the follow- 
ing equation for zq 

l + zoh + zl{l-r)^Q. (23) 

Considering the above as an equation with respect to 
zq its solution could be a complex number on the unit 
circle only if r = (u |(X]u "^i- |0u)(0«|)| w) — 0. As the 
operator R = "^i- |0i;)(0„| = Y.v'^-"^^ positive, 
as well as a sum of positive operators, it follows that 
Pv\u) = |0i,) (0„|u) = for all t; G G to which tails are 
attached. From this it is obvious that the state \u) is a 
bound state of the second kind. Equation reduces 
to 

{{1 + zI)I + zoHg)\u)^Q (24) 
which can be rewritten as 

Hg \u) = - (^zo + \u) . (25) 

Comparing this to ((2T|) we find the energy of the state 
\u) to be Eu = -(zo + 1/zo). 

We want to prove that even when A{zq) is not invert- 
ible as in the case described above, (fTTj) still has a well- 
defined non-zero solution which in addition is orthogonal 
to all bound states. First we prove that if \u) satisfies 
equation then it satisfies A{zq) \u) = as well: 

A{z;) \u) - (/ + zIHg + (z*)2(/ - R)) \u) 
= {l + z*E^ + {zlf)\u) 

= (l - 4 (^0 + ^) + (^0*)') I") = 0. (26) 

The last equality is true because zq lies on the unit circle. 

Let's denote by Kz^ the orthonormal projector on the 
kernel of the operator A{zq) in the Hilbert space of 
G. In other words K^a is the projector on the linear 
subspace spanned by all \u) which satisfy ((22|) . Thus 
has the properties A{zq)Kz„ — and as we saw 
KzgPv — PvKzo — 0. From the above proof we also see 
that A{zl)Kza = 0. Then 

Kz,A{z,) = {A\z,)KlJ = {A{z*,)KzJ - 0. (27) 

From this it follows that / — K^a is the orthonormal 
projector on the image of A{zq) and as (/ — Kzq)Pv' = 
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Pv'{I — Kzg) = Pv' we see that |0t,') is in the image of 
A{zq). Thus a solution to PT|) exists. To define the prop- 
agating state for z = zq we choose the only such solution 
that is orthogonal to kernel of A{zo) which ensures its 
orthogonality to all bound states of the second kind. To 
do this we consider the pseudo-inverse of A(zo) defined 
by 



A-\zo) = {A{zo)+K,„y 



(28) 



which we can use to define a solution to pT|) with the 
necessary properties: 

\zo,vm) = {l^zl)A-^(zo) |0,„) 

= {l^zl){A{z^)+Kz,)-^\%). (29) 

This concludes the proof of the orthogonality of propa- 
gating states and bound states. 



where the sum is over all vertices connected to v. If we 
define = {v\tjj) and take \i'{t)) = \(j){t)) the equation 
takes the more standard form 



dt 



(36) 



The quantity V'i^V'i' is just the probability at vertex v at 
time t and thus (IV') 1 1'0)) can be thought of as the 
probability current at time t fiowing from vertex v to 
vertex v' . 

Let's assume now that |V'(0)) and |</3(0)) are eigenstates 
of with equal energy E. Then 



dt 



= 



and so 



(37) 



(38) 



IV. PROPERTIES OF THE S-MATRIX 



The S-matrix is defined as 



^vm,v' m' — tvm,v' m' 



(30) 
(31) 



are given by ^ . The dimension 
of the S-matrix is equal to the number of tails attached 
to the graph G. To simplify the notation instead of using 
double indices for the S-matrix elements we will use just a 
single one, r, which will stand for the ordered pair (v, to). 
We will also use the notation |r) for the state to which the 
tail labeled by r = {v,m) is connected: |r) = |0t,) = \v). 
From (fTT|) it follows that the elements of the S-matrix 
will be given by 



Srr'{z) = il-Z^){T\A{z)-^\T')-Srr'. 



(32) 



When A{z) is uninvertible, A{z)^^ should be thought of 
as a pseudoinverse in the sense of the previous section. To 
prove the unitarity of the S-matrix we need to introduce 
the following quantity 

fv'AW . \^)) = * {{^\v2){vi\^) - {^\VI){V2\^)) (33) 

where \^p) , \ip) G Hg and the two vertices, vi,V2 G G, 
are connected to each other by an edge. A obvious prop- 
erty of j^j^ that will be used later is the antisymmetry in 
the vertices: = — j"^. Let's consider the Schrodinger 
equation for the graph with the Hamiltonian given by ^ 



dt ~ ' 



G I 



(34) 



Any two of its solutions, \ip{t)) , \(f>{t)), will satisfy the 
following "conservation" equation 



d{'ijj\v) {v\(j)) 
dt 



^j:'(iv'),|0)) = o 



(35) 



for any vertex w G G. Summing the formula above over 
all vertices in G we have 



EE 



(39) 



In the above sum for each term , \ip)) there exist a 

term jl^i{\ip) , \^)) whenever v, v' G G. Because of the an- 
tisymmetry of every such pair of two term will cancel 
out, leaving us with terms for which one of the vertices 
is on a tail: 

E E J'^^'d^) ' 1^)) = E^o;:(iV') , \^)) = 0. (40) 

v£G v' r' 

If we choose — W) — \z,t) (the solution to equation 
^ which when restricted to Tic satisfied equation ([TT])) 
the above equation reduces to 



1. 



(41) 



Here we use the convention the a primed sum signifies 
that we sum over all r' ^ r. If we choose = l^^, ti) 
and \tp) = \z,T2) we obtain 



E' 



0. 



This proves the unitarity of the S-matrix. 
Another property that we will need is 

Stt'{z*) = {s^)Tr'{z) = S*r'r{z)- 

Equation is easy to prove using ([5^ : 

S,r'{z*)^{l-{z*f){T\A{z*)-'\T') 

= (l-z2)*(r|(A(z)t)-i|r') 
= {1 - z^ {t' \A{z)-^\Ty - 5,,, 

= {{I - Z^) {t' \A{z)-^\t) - 5rr') 
= Sr'riz) = (S^)rr'(^)- 



(42) 



(43) 



5tt' 
Stt' 
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V. VERTEX STATE BASIS AND ENERGY 
EIGENVALUE BASIS 

In this section we present formulas for the spectral de- 
composition of the identity corresponding to the Hamil- 
tonian H^. It is easy to obtain after determining its 
eigenstates — the propagating states together with the 
bound states of the first and second kind. For v,v' S G 
we have 

b 13 

+ i /E^'=.-(^)^fc,-(^')'^^' (44) 
- 

where 

ipbiv) = {v\xb), 

In Dirac notation the formula will give the resolution of 
the identity operator I-j-c^ for the Hilbcrt space Hq: 

+ ^/EI^'^)(^'^I^^' (45) 

Using this formula we can easily obtain a simple expres- 
sion for a vertex state |7i^) corresponding to a vertex on 



the tail r of the graph G in terms of the energy eigen- 
states: 

\n,) =Y,{^b\nr)\>ib) ^Y.^(i\nr)\li) 

b 13 

TT 

+ ^ jY.^k,T'\nr)\k,T')dk. (46) 
^' 

This can be further simplified. As we have already 
proven, the bound states of the second kind have zero 
overlap with any vertex that lies on a tail, (/^In,-) = 0. 
Thus the second sum in the above expression is zero. 
Formula (llSp also implies that {Kb\nr) = e~^''"(>fb|Or). 
From (|4|) and ([5]) we see that 

{k,T'\nr) = 5..'e*" + s;,,(fc)e-*'=", (47) 

or in terms of the z variable 

{z,T'\nr)=5rr'Z^ + S%,,{z)z-^, (48) 

Here we need the following identity: 

^s,r{k)\-k,v) = \k,T) , (49) 
or in terms of the z variable 

Y,S.r{z)\z\u)^\z,T). (50) 

First we prove this for the restriction of the propagating 
states on the graph G: 



^ \z\vf = E(l - ^*^)A-Hz*) W) ((1 - z^) (v \A-\z)\t) - <5.0 

= ((1 - ^')RA-\z) - I) \r) 

= (1 - z^)A-^\z) {z*^A{z) - z*2(l z^)R) A-^{z) |r) 

(1 - z'^)A-^\z) + z*H + I-R- z*^R + R) A-^{z) \t) 
= (1 - z^)A-^\z)A^{z)A-\z) |r) = |z,t)^ (51) 



In the above proof we used that \z\ = 1. All inverses 
above are well defined as well because they act on |t) 
which lies in the image of A{z). 

Now we need to prove the identity for any vertex lying 
on a tail. Multiplying ([50)1 by jn^./) we arrive at the 
following formula which we need to prove: 

^s^riz){nr'\z*,iy) = {nr'\z,T). 



It is easily proven using and ([^5]) : 

^S^r{z){nr'\z*,iy) = ^S^r{z){z*,l'\nr')* 
= S„r{z){Sr't.zJ' + Sr'^{z*)z-'') 

= ^ s,,(z)(<5,,,z" + (st),,,(z)z-") 

= Sr'r{z)7j' + (5.'^^"" = (<5r'r2" + 4,,(z)z-")* 

= {z,T\nr>)* = {nr'\z,T). (52) 
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This proves (|49)) . 

Now we ean simplify the sum under the integral in 
(gni) with the help of (gT]),!!!]) and the unitarity of the 
S-matrix: 



^(fc,rV,)|fc,T') 

r' 

r' 

= e*'^'"|/c,r)+e-**="|-/c,r). 
Substituting this back into pS]) leads to: 

b 

+ ^ j{e'^^%T)+e-^^-\-k,T))dk 



= ^(z,K>|z,) + i- /■ e^^-\k 

I. 'J 



dk. 



(53) 



Using the z variable this formula takes very appealing 
form: 



^^)=J2z^{z,\0^)\z,) + ±-<|^z-\z,r)^. (54) 



VI. CUTTING A TAIL 

In this section we want to investigate the effect of cut- 
ting a tail on the S-matrix of the graph. We will show 
that the S'-matrix of the new, pruned graph can be ex- 
pressed in terms of the elements of the S'-matrix of the 
old graph. Thus we are given a graph G to which n tails 
are attached and we assume that we know its S-matrix. 
Without loss of generality we can choose to cut the first 
tail which we denote by Tc- For the unpruned graph and 
pruned graph the elements of the S-matrix are given by 
([5^ with 



A{z)=I + zHg + z^Q, 
Aciz) = I + zHg + z^Qc 
= A + z^\tc){tc\, 



A + z'Pr, 



(55) 
(56) 



respectively. As can be seen from ([5^ the problem re- 
duces to expressing the matrix elements of A~^ in terms 



of the matrix elements oi A ^ . 

{r\A-^\T') 

= {T\{A + Z^PrX'W) 

= ij\[I + 7?A-^PrX^A-^\T') 

oo 

= {T\Y^{-7?A-^p^yA-^\T') 

OO 

x^(-z2 (r,|A-i|r.))^-i(r,|A-i|r') 
i=i 

"^"'^ '"^ l + z2(r,|A-i|r,) ■ 

Multiplying both sides of this expression by 1 — and 
using ([5^ we get 



1 + Z Sr^Tc 



(57) 



A. Leaving a Stump 

In stead of cutting the tail at the root now we want 
to leave a stump of length L. This case can be easily sub- 
sumed in the previous one. We define a new tail fc which 
coincides with Tc but its beginning is at the L-th vertex of 
the tail Tc, \Lt^- Thus |0f^) = 1^,-^) The propagating so- 
lution to ([2) in this case will be equal up to a phase to the 
solution when we think of the tail being attached to the 
root |OtJ- The conditions that the reflection and trans- 
mission coefficients need to satisfy is given by ^ and 
([5]) which is only possible when the phase between the 
two solutions is chosen appropriately: \k)~ — e^^^ \k)^ . 
Thus 



ikn 



= e 
(rv|fc}f, 



z ~ij -J- ^ ^ • 



For energy eigenstates with the same energy but which 
are incoming on a different tail we get 



{nfjk)r = {{n + L)rjk)r 



z^U 



All other such conditions that don't involve the tail that 
is being cut are satisfied automatically. From the above 
equations and the definitions of the reflection and trans- 
mission coefficients (H)) we see that the elements of the 
S-matrix for the graph with tail fc satisfy: 



Z^^S 

•2 ^TTc -i 
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Plugging those in (|57|) we finally get: 



Ta and f. After some simplification we get 



,2(L+1), 



l + z2(i+i)s^^^ 



(58) 



B. Cutting k Tails 

It is obvious that if we cut more than one tail it doesn't 
matter the order in which we cut them. Let us assume 
that the 5-matrix has the following block form: 



Vkx(rn-k) 



Wkxk 



(59) 



(1 - z^)s^r, 

(1 ~ z'^)Sr^r 

(l-z2)(l + s.,^,) 
X "4" s -7- ~ -J- - 



(61) 



Here m, as defined in Chapter [III is the total number 
of tails and k is the number of tails that will be cut. 
Thus we want to cut the tails corresponding to the last fc 
entries in the S'-matrix. We obtain the following formula 
for the new S'-matrix: 

S'' = T-z^U(l + z'^W)-^V. (60) 

(The identity matrix in the above formula is a fc x fc 
matrix.) It is obvious that formula ([57]) follows from the 
above. 



VII. ATTACHING A TAIL 



VIII. CONNECTING TWO TAILS TO FORM 
AN EDGE 

The setup is as in the previous sections but the goal 
now is to connect two tails — in other words to cut two of 
the tails connected to the graph and replace them with an 
edge between the vertices they were connected to. Again 
without loss of generality we connect the first and second 
tails, Ti and ti. When we do that the Hamiltonian and 
the operator Q change to 

Hg = He — X(^2)i 

Q = Q + P(2) 



Let's assume that we want to attach one more tail to 
a vertex that already has at least one tail attached to 
it. We want to express the S'-matrix of the new graph 
in term of the old one. Of course, it is also possible to 
add a tail to a vertex that doesn't already have one, but 
in this case the new S'-matrix cannot be derived simply 
from the S'-matrix of the original graph. Let us denote 
the new tail being attached with t^, the vertex we attach 
it to with V and the tail already attached to the vertex 
V with T5. We note that Ir^) = \Ta) = \v). By similar 
argument to the one we made in section IVII we get 



{t\A-^\t') = {t\A-^\t') 



(t|A-1|t,) {rAA-^\T') 



where 



An — A Z Pt, 



A 



where we have used the notation 

X(2) = |n)(r2| + |r2)(Ti|, 
P(2) = |n)(ri| + |r2)(T2|. 

This leads to the following expression for A{z) 



i = / + zHg +Z^Q^ A- ZX(2) + Z^P^2). 



For convenience we denote B = z'^P(2) ~ ^^(2)- Again 
we look at the matrix elements of A^^: 



(T|i-i|r'} - {t\{A + B)-^\t') 

00 

= (r| (/ + A-^B)-^A-^ \t') = (r| ^{-A'^ByA-^ |r') 



In order to reexpress everything in terms of the elements From the definition of B it follows that 
of the old S'-matrix we consider different cases depending 
on the first and second indices of the S-matrix. In the 

below r and t' signify tails that are different from both -^(2)^-^(2) = P{2)B = BP(2) = B, (62) 
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from which we see that 



{r\A- 

OQ 

= 5](r|(-A-%)BP(2)FA-Mr') 



(r| A-iP(2) J2 B{-P^,)A-'P^2)By-'P^2)A-' |t') 



- (t| A-iP(2)B(P(2) + P(2)A-iP(2)P)-iP(2) |r') 

- (t| A-iP(2)(P-i + P(2)A-iP(2))-iP(2)A-i It') . 

(63) 



In the above whenever necessary the inverses should be 
thought of as pseudo-inverses. 

We want to use a block-matrix representation again. 
We will use an almost identical notation as the one we 
used before 



5* = 



Tn X n Un x 2 
^2xn W^2x2 



(64) 



but now the last two entries in the S'-matrix correspond 
to the tails that are to be connected. Multiplying 



by 1 — z and after some simplifications we get: 
S - zU{zW - X)-^V 



(65) 



with X being just the X Pauli matrix. 

To give an explicit formula for the elements of the 
new S'-matrix we need the pseudo-inverses of B^^ + 
P(2)A-^P(2): 



B- 



D{z) 



(66) 



where 



D{z) 

C(2) 



: 1 — (ai2 + a2i)z + (an + 022)2^ 

- (aiia22 - 012021)2^(1 - z'^), 

- 022 |ti)(ti| - ai2 |Ti)(r2 
- 021 1^2) (ri I + an \t2){t2 \ , 

where defined through the following formula: 

P(2)A"^P(2) = an \ti){ti\ + ai2 \ti){t2\ 

+ 021 |t2)(ti| -f 022 |r2)(r2| . (67) 



Substituting in ([55)1 and again using the definition 
of the elements of the S-matrix ([5^ we find 



z(^S'rTi^r2T' ~t- 5-J--J-2 S-T-^ 7-' ) 



A. Composition of unitary gates 

The universality of quantum walks in continuous time 
was proved in In the model presented there a quan- 
tum wire corresponds to a set of tails connected to a 
graph. Specifically the state of a qudit is represented by 
a linear superposition of incoming waves on d of the tails 
with the same energy. Thus to each tail corresponds one 
of the orthogonal states of the qudit. The graph to which 
the tails are connected will implement the quantum gate. 
Another set of d tails connected to the graph represent 
the quantum wire which carries the state of the qudit 
with the quantum gate applied to it. The order for this 
to represent a quantum gate it is needed that a wave com- 
ing in on any incoming tail needs to scatter on outgoing 
tails only for the fixed energy at which the computation 
is performed (the scattering matrix of a graph cannot 
satisfy this condition for every energy). Thus it is easy 
to see that for that particular energy the S-matrix for 
this graph needs to have a block form: 



S,o 

So^ 



(69) 



where the first d entries stand for incoming tails and the 
last d for outgoing. Thus the quantum gate being imple- 
mented is given by the unitary matrix Soi- 

We want to show that if we connect two graph using 
the rules for connecting tails we will end up with a gate 
that is composition of the quantum gates corresponding 
to each graph. Before we have connected the two graphs 
the S-matrix is just the direct product of the S-matrices: 



5" = 



/ Si 

Sl^ 

Si 

V si 



(70) 



We will use a generalization of formula (|65p to find the 
S-matrix after the connections are being made. The for- 
mula retains the same form even if we connect more than 
two tails as long as we arrange the tails being connected 
to correspond to entries in the lower right corner of the 
S-matrix. The generalization X of X in formula (|65|1 is 
going to be given by either 



X = 











M 





X 



















V 











(71) 



-I- 



1 - 




^Ti T2 '^T2 Tl 


)z2 


or 




1 - 


K'^rir2 ~t~ '^T2Ti)^ (-^nri 'Sr2T2 ^ 


r2 ^T2Ti 


)z2 




X 



/ 
/ 



(72) 



^TiT2 ■^T2Ti 



depending on how the tails being connected are arranged 

(68) 

in the matrix W . Thus by permuting the entries (I70p we 
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obtain: 



lemma. 



S' 




sl\ 
si. 


/ 



(73) 



Thus for the block matrices T, U, V and W we have 

T = 0, C/=(4 f ), (74) 
^ = ( Jl. ^° ) , W^O (75) 



and X needs to be of the second form. For the new S- 
matrix we obtain: 



S - zU{zW ~ X)-^V 



s}^ 
si 







-/ 



sf^ 
si 





:2 oi 



'-'io'-'io 



'-'oi'-'oi ^ 



(76) 
(77) 

(78) 



We see that that unitary being implemented is given by 
^SoiSoi which is exactly equal to the composition of the 
two unitaries up to the phase z. The edges that connect 
the two graphs which are formed after connecting the 
tails lead to the introductions of this phase. 



IX. UNITARITY PRESERVATION IN CUTTING 
AND CONNECTING TAILS 

In this section we prove that the above operations pre- 
serve the unitary of the S-matrix. We need the following 



Lemma. IX. 1 Given a unitary matrix S with block- 
form representation 



S 



Vm X n Wm x m 



and an unitary matrix Gmxm such that G + W is invert- 
ible, the matrix 



s -u{G + wy^v 



(79) 



is unitary. 



Proof. The unitarity of S implies that 



TV^ = ~UW\ 
VV^ = I -ww\ 



T^T = 1 - y V, 

VT^ = ~WU\ 
U^U = I -W^'W. 



(80) 



In the above formulas, although not indicated, the iden- 
tity / needs to be understood as having the appropriate 
dimension for each equation. Then 



ss^ = {T- u{G + wyW) (T - u{G + wyWy 

= TT^ - U{G + Wy^VT"^ - TVHG^ + W^y^U'^ + U{G + Wy^VV^{G^ + iyt)-i[/t 

= / - uu^ + u{G + wy^wu"^ + uw\g^ + w^y^u^ + u{g + wy^{i - ww'^){g^ + w'^y'^u'^ 

= i-u{i-{G + wy^w - iyt(G't -f- M/t)-i - (G -I- wy^iG^ + w'^y^ + (g + wy^ww^G'^ 4- w'^y^) t/^ 

= / - c/ ((/ - (G + wy^w){i ~ w\G^ + + wy^iG'^ + VKt)-i) 

= I-U{{G + wy^GG^G^ + W^y^ - (G + wy\G^ + Tyt)-i) [/t = J (81) 

I 



This proves the lemma. 

In the case of cutting a tail from formula ((60)) we see 
that G = I/z^ which is unitary because \z\ — 1. 

In the case of connecting two tails from formula ([65]) 
we see that G = —X/z which is obviously unitary as well. 

This proves the unitarity of the S'-matrix after an op- 
eration of cutting or connecting tails. 



X. DISCUSSION 

Quantum walks in their many varieties have proven 
immensely useful in the area of Quantum Computation 
as a tool for developing new algorithms or as an abstract 
model to studying the behavior of quantum systems. The 
model we are concerned with in this paper involves quan- 
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turn walks on infinite graphs, it is natural to develop 
scattering theory for these systems. We have studied 
both propagating and bound states, and how they de- 
pend on the structure of the finite graph to which the 
tails are attached. The equations that define the propa- 
gating and bound states are an example of the quadratic 
eigenvalue problem. The S-matrix, its unitarity, defini- 
tion and conservation of probability density and current 
were addressed, as well as the orthogonality of propagat- 
ing and bound states. We were able to derive formulas 
for the S-matrix of a graph that is obtained by opera- 
tions on the tails in terms of the S-matrix of the original 
graph. 

In future work, we would like to pursue the defini- 
tion of hitting time for such graphs. In previous work, 
0, m, [13, [3 : properties of the hitting time for finite 
graphs were explored. We would like to give appropriate 
definition for hitting time when we consider the tailed 
graphs. One such definition may follow in the footsteps 
of [l8| , where the quantum walk starts from an arbitrary 
state, and we are interested in the probability and av- 
erage time to hit any of the tails. A different approach 
could be followed if we consider a wave packet that is in- 
coming along one of the tails with an energy in a narrow 
band. Then a definition for a hitting time may be the 
time it takes for the packet to reach a different outgoing 
tail to some standard graph. For example, if we have a 
graph with just two tails, we can compare the time for 
which a packet travels through the graph from one of the 
tails to the other and compare it with the time that the 



same packet will take to move along the infinite line. 

Defining the hitting time will permit us to explore new 
ideas and paradigms for quantum algorithms. In [l3] in- 
finite hitting times were discovered for the discrete-time 
case, and later the same phenomenon was observed for 
the continuous-time case in |il8]. The existence of an 
infinite hitting times can lead to possible quantum al- 
gorithms. The algorithms can be based on the simple 
observation that the quantum walk may or may not hit 
a certain set of vertices depending on the existence of in- 
finite hitting time. This would serve as the final result in 
our quantum computation, much like the evaluation of 
the NAND tree is based on whether the quantum walk 
hits the second tail or not. Both infinite hitting times 
and exponentially fast quantum algorithms as in [1^ are 
based on the idea that the quantum walk is constrained 
to evolve on a subspace of the whole Hilbert space be- 
cause some special property of the graph. Symmetry of 
the graph is one property that can lead to these effects, 
but other properties may do so as well. This is an open 
problem, and an area of active research. 
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